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ABSTRACT 

We present extensive radio observations of SN 2003L, the most luminous and 
energetic Type Ibc radio supernova with the exception of SN 1998bw (associated 
with GRB 980425). Observations from the Very Large Array are well described 
by a fitting a synchrotron self-absorption model to the emission spectrum. This 
model implies a sub-relativistic ejecta velocity, v ~ 0.2c, and a size of r pa 
4.3 x 10 15 cm at t ~ 10 days. The circumstellar density is suitably fit with a stellar 
wind profile, n e oc r~ 2 and a constant mass loss rate of M pa 7.5 x 10 -6 M Q yr _1 . 
Moreover, the magnetic field follows B oc r _1 and the kinetic energy of the 
radio bright ejecta is roughly E pa 10 48 erg assuming equipartition of energy 
between relativistic electrons and magnetic fields. Furthermore, we show that 
free-free absorption does not contribute significantly to the radio spectrum, since 
it implies ejecta velocities which are inconsistent with size constraints derived 
from Very Long Baseline Array observations. In conclusion, we find that although 
SN 2003L has a radio luminosity comparable to that seen in SN 1998bw, it shows 
no evidence for a significant amount of energy coupled to relativistic ejecta. Using 
SN2003L as an example, we comment briefly on the coupling of ejecta velocity 
and energy in Type Ibc supernovae. 

Subject headings: gamma rays: bursts - radiation mechanisms: nonthermal - 
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1. Introduction 

Radio emission from core-collapse supernovae (SNe) have been sought primarily to study 
the circumstellar medium (Weiler et al. 1986; Chevalier 1998). However, starting with Type 
Ibc SN 1998bw, associated with the under-luminous gamma-ray burst GRB 980425 (Pian et 
al. 2000; Galama et al. 1998; Kulkarni et al. 1998), radio observations have also been used a 
tracer of relativistic ejecta (Kulkarni et al. 1998; Li & Chevalier 1999). This event prompted 
a storm of publications: SN 1998bw was the most luminous radio supernova (at early times), 
showed strong variability and produced copious gamma-ray emission relative to ordinary 
SNe Ibc. These data provided the first observational evidence that at least some SNe Ibc 
are associated with GRBs (e.g. Bloom et al. 1999, 2002; Stanek et al. 2003, c.f. Price et al. 
2003). 

To further explore the connection between GRBs and SNe Ibc, we began a radio survey 
of local Type Ibc supernovae. Our survey was motivated by two simple phenomenological 
questions: (1) What is the prevalence of SN 1998bw-like supernovae?, and (2) Is there a 
continuum between ordinary SNe and GRBs as traced by their relativistic ejecta? 

In Berger et al. (2003) we reported the first three years of the survey. We found that 
most SNe Ibc are not detectable with current radio sensitivity, with limits reaching 10~ 3 
that of the peak brightness of SN 1998bw. Clearly, luminous SN 1998bw-like supernovae are 
rare (less than 3%) and those within an order of magnitude of SN 1998bw are uncommon 
(less than 10%). 

Here, we present radio observations of SN2003L, the first SN Ibc with a radio luminosity 
comparable to SN 1998bw. While there is no evidence for mildly-relativistic ejecta, the radio 
emitting ejecta of SN2003L are unusually energetic. This discovery suggests the existence of 
a sub-class of SNe Ibc with similar energetics to SN 1998bw and even weaker central engines, 
if any. 

The organization of this paper is as follows: observations from the Very Large Array 
(VLA), the Very Long Baseline Array (VLB A) and the Chandra X-ray Observatory (CXO) 
are described in §2 and §3. Preliminary estimates of the energy and velocity are presented in 
§4. Modeling of the radio light-curves is discussed in Sections 5 and 6 assuming the dominant 
absorption is from internal synchrotron self-absorption and external free-free absorption, 
respectively, Modeling of the X-ray emission follows as §7. In the Discussion (§8) we conclude 
by addressing the nature of SN 2003L. 
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2. Radio Observations of SN 2003L 

2.1. Very Large Array Data 

SN2003L was optically discovered on 2003 January 12.15 UT, offset 9".0 W and 1".5 
N from the center of NGC 3506 (Boles 2003) and at a distance of d « 92 Mpc (Valenti et 
al. 2003). In our first observation with the Very Large Array 7 (VLA), on 2003 January 
26.2 UT, we detected a radio source coincident with the optical position at a(J2000) = 
ll h 03 m 12.3 s ,<5(J2000) = + 11°04'38.1" (±0.1 arcsec in each coordinate) with flux density of 
/„ = 743 ± 39 /iJy at 8.5 GHz. This, and subsequent observations at 4.9, 8.5, 15.0 and 22.5 
GHz, are summarized in Table 1. 

All observations were taken in standard continuum observing mode with a bandwidth of 
2 x 50 MHz. We used 3C286 for flux calibration, and for phase referencing we used calibrators 
J1118+125, J1120+134 and J1103+1 19. Data were reduced using standard packages within 
the Astronomical Image Processing System (AIPS). 

As seen in Table 1, observations were taken in each of the four VLA configurations. The 
decreased spatial resolution of the C and D configurations resulted in detection of diffuse 
host galaxy emission surrounding the supernova, primarily affecting the low (< 8.5 GHz) 
frequency observations (Figure 1). At 8.5 GHz the host galaxy extends 45 x 45 arcsec, but 
is resolved away on baselines > 10 kA. By restricting the UV range to exclude shorter 
baselines, we measured the flux density of the supernova with minimal contribution from 
the host galaxy. Measurements at higher frequencies and in other VLA configurations were 
made by fitting a Gaussian to the SN emission and solving for the integrated intensity. 

Figure 2 shows the SN 2003L radio light-curves assuming an approximate explosion date 
of 2003 January 1.0 UT as derived from optical data (Soderberg et al. 2004). The radio 
light-curves evolve as f u oc t 1 - 2 before the peak, exhibit a broad maximum and fade as 
/„ oc t~ L2 . The spectral index ((3 with f v oc u 13 ) is observed as (3 « —1.1 in the optically 
thin regime (Figure 3). 



7 The Very Large Array and Very Long Baseline Array are operated by the National Radio Astronomy Ob- 
servatory, a facility of the National Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 
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2.2. Very Long Baseline Array Data 

Complimentary to our VLA campaign, we obtained a single observation of SN2003L 
using the Very Large Baseline Array (VLBA). The six hour observation began on 2003 March 
7.30 UT (t ~ 65 days), and was taken in standard continuum mode with a bandwidth of 
4x8 MHz centered on the observing frequency of 4.9 GHz. Fringe calibrations were applied 
using 3C273 and phase referencing was conducted using J1103+1158 at an angular distance 
of 0.9 degrees from the SN. 

We detect SN2003L at a position of a = ll h 03 m 12 s .3008, 5 = +11°04'38".08 (J2000.0) 
with a positional uncertainty of 10 mas in each coordinate (Figure 4). We note that these 
errors are dominated by the positional uncertainty of J1103+1158, assumed to be at a = 
ll h 03 m 03 s .5299, 5 = +11°58 , 16 /, .61 (J2000.0). Using the VLBA utilities within AIPS, we 
find a flux density for the supernova of F4.9 gh z = 848 ± 64 fiJj within the 3.158 x 1.285 
mas beam. We note that there is no emission from the host galaxy at this resolution. At a 
distance of 92 Mpc, this unresolved VLBA detection places a direct constraint on the size 
of the expanding ejecta of r < 9.1 x 10 17 cm and an average expansion speed of v < 5.4c 
during the first 65 days after the explosion. 

3. X-ray Observations with Chandra 

We observed SN 2003L with the Chandra ACIS-S detector beginning at 2003 February 
10.15 UT (t ~ 40 days). During the 30 ksec observation, we detected 30 ± 6 counts from 
the supernova (Kulkarni & Fox 2003). For spectral extraction we adopt a source aperture of 
2.46 arcsec, which encloses approximately 98% of the flux of an on-axis point-source at 1 keV 
(the peak of our raw photon counts spectrum). Our aperture correction is therefore minimal 
by comparison to our ~20% flux uncertainties (see below). Background counts are selected 
from an annular region extending from 2.46 arcsec to 7.38 arcsec from the source; in general 
larger background regions are preferred for Chandra data analysis but the extended emission 
from the host galaxy makes a more local selection preferable in this case (see Figure 5). 

Schlegel, Finkbeiner & Davis (1998) dust maps give E(B — V) — 0.021 mag for the 
location of SN2003L. Dickey & Lockman (1990) suggest an average N H of 2.25 x 10 20 cm- 2 , 
while a direct conversion from the Schlegel et al. optical reddening suggests 1.13 x 10 20 cm~ 2 
(Predehl & Schmitt 1995). We fix N u for our fits at 2.5 x 10 20 cm- 2 . 

For a power-law spectral fit we find the power-law photon index T = 1.53 ± 0.4 (range 
is 0.86 to 2.24 at 90% confidence). For a thermal bremsstrahlung spectral fit we find the 
plasma temperature has kT > 1.7keV at 90% confidence, with a best-fit value of 6.7keV 
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that is unconstrained from higher energies. The relatively small number of counts enables 
satisfactory fits (reduced chi-squared values of xl ~ 1-3; 5 degrees of freedom) for both 
models, with the neutral hydrogen absorption applied to either a power-law or thermal 
bremsstrahlung continuum. 

Our source count rate is 1.0 ± 0.2cts/ksec in the ACIS-S3 detector. This corresponds 
to a 0.5-5.0 keV flux for our spectral models of approximately 5.3 x 10~ 15 erg/cm 2 /s in 
either case. Extrapolating to the full 2-10 keV X-ray band, the corresponding fluxes are 
9.2 x 10~ 15 erg/cm 2 /s for the power law model and 7.6 x 10~ 15 erg/cm 2 /s for the thermal 
bremsstrahlung model. The associated X-ray luminosity is thus L2-10 kcv ~ 9.2 x 10 39 erg/s 
and 7.2 x 10 39 erg/s for the power-law and thermal models, respectively. In comparison with 
other X-ray supernovae, SN 2003L is among the most luminous ever detected, only a factor 
of ~ 10 fainter than SN 1998bw at a comparable epoch (Pian et al. 2000). 



We are interested in the energy and expansion velocity of the ejecta producing the radio 
emission. As a preliminary constraint, we estimate the brightness temperature, Tg, of SN 
2003L and compare it with robust constraints imposed by the Inverse Compton Catastrophe 
(ICC; Kellermann & Pauliny-Toth 1981). The brightness temperature of a source with 
angular radius, 9 is given by 



Compact sources with T B > 10 12 K cool rapidly via inverse Compton scattering. As an initial 
estimate for the physical size of the supernova, we first assume that the optical expansion 
velocity of 12000 km s _1 (Matheson et al. 2003) can be used as an average speed to describe 
the motion of the radio bright ejecta. Using our approximate explosion date, we estimate the 
shock radius to be r pa 2.9 x 10 15 cm at t ~ 28 days. For an observed flux density of f u pa 3.1 
mJy at v pa 22.5 GHz, the brightness temperature is Tg ~ 6.3 x 10 11 K, just below of the 
inverse Compton catastrophe (ICC) limit. This exercise provides no evidence for relativistic 
ejecta, therefore suggesting that SN 2003L expands with a modest sub-relativistic velocity. 

For radio sources dominated by synchrotron-self absorption (SSA), the brightness tem- 
perature can be further constrained to Tg < 4 x 10 10 K. By assuming equipartition between 
the energy in electrons (e e ) and magnetic fields (e#), we derive an estimate of the velocity 
and a lower limit on the kinetic energy of radio emitting material. This approach has been 
used extensively in discussion of extra-galactic sources and was more recently applied to su- 
pernovae (Chevalier 1998; Kulkarni et al. 1998). The method requires three observables at 



4. Preliminary Constraints 




(1) 
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any single epoch: the spectral peak frequency, is p , the peak flux density, f p , and the spectral 
index of the optically thin emission, f u oc u 13 . Here, (5 = — (p — l)/2 for electrons accelerated 
into a power-law distribution, Nf^y) oc 7 _p above a minimum Lorentz factor, r ) m . Following 
Kulkarni et al. (1998) and Berger, Kulkarni & Chevalier (2002), the angular radius of the 
radiosphere, 9 ep , and the equipartition energy, E cp , are given by 




At t ~ 30 days (our first epoch) the peak flux density is f p ~ 3.2 mJy at v v w 22.5 
GHz. Using the observed optically thin spectral index of /3 ~ —1.1, Equations 2 and 3 
give 8 ep ~ 7.7 yuas implying an equipartition radius of r cp ps 1.1 x 10 16 cm and an average 
expansion speed of v cp ~ 0.15c. The minimum energy of the ejecta is thus E ep w 1.5 x 10 47 
erg and the magnetic field is B cp ps 1.1 G. These constraints alone demand that SN2003L 
is among the most energetic Type Ibc supernovae in terms of radio-emitting ejecta, second 
only to SN 1998bw. We note that additional absorption terms (e.g. free-free absorption) and 
the inclusion of shocked protons serve to increase the equipartition energy estimate. 

Equipartition analyses may also be used to roughly constrain the minimum electron 
Lorentz factor and the characteristic synchrotron frequency, v m . Equating particle kinetic 
energy across the shock discontinuity, we find 



'v\ 2 

rripC 2 7m e c 2 (4) 



where 7 is the average Lorentz factor of the electrons and v is the velocity of the ejecta. 
Assuming roughly 7 ps 7 m and substituting Equation A6 for j m , we derive the following 
estimate for the characteristic synchrotron frequency, u m , 



Adopting the values derived from equipartition analyses (v ps 0.15c, B = 1.1 G and e e = 0.5), 
we find roughly v m ps 1 GHz at t ~ 30 days which is below our observing band. 
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5. Internal Synchrotron Self- Absorption 

In the Appendix we present a rigorous formulation of the temporal and spectral evolution 
of synchrotron emission arising from sub-relativistic supernova ejecta. This prescription, 
based on the formalism of Frail, Waxman & Kulkarni (2000), is generalized to include the 
cases where v m , is greater than the self- absorption frequency, v a , and where e e ^ e B . 

The observed flux density and radio spectrum at any single epoch are determined by 
three parameters: Cf, C T and z/ m ,o- Here, Cf and C T are normalization constants of the 
flux density and optical depth, respectively, while u m is the value of u m at epoch t . The 
parameters Cf, C T and z/ m are in turn determined by the values of four physical parameters 
at t = t : the magnetic field, B , the shock radius, r , the minimum electron Lorentz factor, 
7m,o, and the ratio of energy densities, do = £ e /tB- With four physical parameters (B , r , 
7m,o, do) and only three constraints (Cf, C T , v m ,o)i we must assume an additional constraint 
in order to find a unique solution. An additional constraint is obtained by adopting a 
value for do- Inverting the equations for Cf, C T and z/ m0 (Equations A13, A14 and A15, 
respectively) we derive the following expressions for B , r and 7 m ,o- 



B 



19^21 ^4 "1 V 1 ? 



1.3 x10V 5 (2 + p)- 6 G D -2)- 4meC V 



r = 



e 17 dod 4 

1.5x10- 2 tt- 9 (2 + p) 7 (p-2)- 1 4^ 



x 



7m,0 



3.1x10-', C + rt'fr- 2 *' 8 ^ 



1/17 



Cj 2/17 C 4 T /17 u m 2 t 2)/17 G (6) 

(7) 



x Cf 17 cr^t 2)/U cm 



x 



r ,l/17 r ,-2/17 I/ (13+2p)/34 



(8) 



Here, we define d as the distance to the supernova and r\ characterizes the the thickness of 
the radiating electron shell as r/rj. We make the standard assumption of a thin shell with 
7] = 10 (Li & Chevalier 1999; Frail, Waxman k Kulkarni 2000). 



5.1. Hydrodynamical Evolution of the Ejecta 

Evolutionary models governing the hydrodynamics of the SN ejecta provide additional 
constraints. As discussed by Chevalier (1996), there are several different models that can be 
used to describe the hydrodynamic evolution of the ejecta. These models allow the temporal 
indices, a r , as, « 7 and a$ to be constrained. Here, we adopt the standard model of Chevalier 
1996 for our basic SSA fit. 
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The first assumption of the standard model is that the hydrodynamic evolution of 
the ejecta is self-similar across the shock discontinuity. This implies r oc t ar with a r = 
(n — 3)/(n — s) where n is the density profile of the outer SN ejecta (p oc r~ n ) and s is the 
density profile of the radiating electrons within the shocked circumstellar medium (n e oc r~ s ). 
In addition, the standard model assumes that the magnetic energy density (Ub oc B 2 ) and 
the relativistic electron energy density (U e oc n e ^ m ) scale as the total post-shock energy 
density (U oc n e v 2 ). In this scenario, the magnetic field is amplified by turbulence near the 
shock discontinuity, implying fixed energy fractions, e e and e#, and thus a constant value of 
5 throughout the evolution of the ejecta. 

Adopting this evolutionary model, the temporal indices are constrained as follows. The 
self-similar solution requires that the total energy density scales as U oc t an£+2av where a ne 
defines the temporal evolution of the electron density. For a radial dependence of s, we 
have a He = —sa r . Using the scaling for the shock velocity (Equation A2), we then have 
U oc £- SQ ?-+2(a r -i) _ Similarly, U e oc t~ sar+a ~* . In a model where e e and e# are constant, it 
follows that a$ = and U e oc Ub oc U . This results in the following constraints, 



a 7 = 2(a r - 1) (9) 
(2 - s) 

«s = L y i «r-l. (10) 

With the standard assumption that the circumstellar medium is characterized by a wind 
density profile, s — 2, Equation 10 simplifies to as = — 1. 



5.2. SSA Model Fit for SN 2003L 

Using the SN2003L multi-frequency radio light-curves and a chosen value of we fit 
for the constants C/ and C T as well as the temporal index a r . As discussed in Section 4, u m 
is estimated to be below our radio observing band and we therefore adopt v m $ 1 GHz. We 
note that the uncertainty in z/ m corresponds to minimal uncertainty in the derived physical 
parameters as shown by Equations 6-8. 

Adopting a reference time of t = 10 days, we find a best-fit solution (xl — 7-5; 105 
degrees of freedom) for parameters: Cf = 7.2 x 10~ 53 , C T = 4.5 x 10 38 , and a r = 0.96. 
According to Equations 9 - 10, these values imply a 7 = 0.075 and a B = —1.0. Here, we 
have used 5o — 1 with a$ = and an electron energy index, p = (—2(3 + 1) = 3.2, based 
on the observed optically thin spectral index, (5 ~ —1.1 (see Section 2 and Figure 3). As 
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discussed in the Appendix, we parameterize the sharpness of the v a spectral break with 
C = [0, 1] and find ( ~ 0.5 for our best-fit solution. The model provides a reasonable fit to 
the light-curves as seen in Figure 2. We note that the large xl is dominated by scintillation, 
primarily affecting the lower frequency observations. 

Better fits can be obtained by relaxing the assumptions of the standard model (hereafter 
Model 1) as described in Section 5.1. In Model 2, we remove the assumption of a wind 
density profile and effectively fit for s. For this model we find a best-fit solution (Xr = 5-7; 
104 degrees of freedom) for C } = 1.2 x 10~ 52 , C T = 2.5 x 10 38 , a r = 0.85 and a B = -0.84. 
By Equation 10, this solution implies a shallow CSM density profile of s = 1.6. As in the 
case of Model 1, we have adopted Jo — 1) ®$ — 0, p = 3.2 and ( = 0.5 for this fit. 

A similar fit is obtained by fixing the density profile as s = 2 and removing the assump- 
tion of constant energy density fractions, e e and eg. In this case, Model 3, the constraints 
given by Equations 9 and 10 do not apply since a$ is non-zero. By definition, # oc U e /UB, 
resulting in the constraint 

a$ = — sa r + a 7 — 2as (11) 

which gives a$ = —2a r + a 7 — 2a b for s = 2. For Model 3, we find a best-fit solution 
(X 2 r = 5.7; 103 degrees of freedom) for C f = 1.2 x 10~ 52 , C T = 2.4 x 10 38 , a r = 0.78, 
a B = —1.1 and a 7 = 0.16. The implied evolution of J is given by a$ = 0.88 with e e oc t 61 
and €b oc t~ - 27 . Consistent with Models 1 and 2, # — 1, P — 3.2 and ( = 0.5 were used for 
this fit. 

A comparison of the best-fit solutions from Models 1-3 are shown in Figure 2. It 
should be noted that Models 2 and 3 produce the same light-curve fit because since they 
both allow for one extra degree of freedom. While these two models provide a slightly better 
representation of the SN 2003L data, we adopt Model 1 as our best SSA solution since these 
assumptions have been shown to be consistent with the full physical model for SN 1993 J 
(Fransson & Bjornsson 1998), currently the best-studied radio SN. 

5.3. Physical Parameters 

Using our solution for Model 1, we find the following values and evolution for B, r and 
7 m . The temporal evolution of the magnetic field is B 4.5(t/to)~ LO G with an associated 
radial dependence, B oc r~ L04 , similar to that found for the Type lib SN 1993 J (B oc r _1 ; 
Fransson & Bjornsson 1998). The shock radius is well described by r ps 4.3 x 10 15 (t/t ) ' 96 cm 
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implying a sub-relativistic velocity of v pa 0.20(t/t )~ ' 04 c. We note that this radial evolution 
is faster to that of SN1983N (r oc t ' 86 ; Chevalier 1998) and is within the expected range, 
0.67 < a r < 1.0 (Chevalier 1996, 1998). The minimum Lorentz factor of the electrons follows 
7 m pa 8.9(t/t )~ ' 08 and the implied evolution of v m (see Equation A10) is v m pa 1.0(i/i ) _1 ' 16 
GHz which is nearly comparable to the evolution observed for SN1998bw (y m oc t~ 10 ; 
Li & Chevalier 1999) and significantly slower than the decay rate in the adiabatic Sedov- 
Neumann- Taylor solution (v m oc t~ 3 ; Zeldovich & Raizer 2002; Frail, Waxman & Kulkarni 
2000). 

Using these derived values, additional physical parameters can be computed. The num- 
ber density of emitting electrons within the circumstellar ejecta, n e , is thus 



n P = 



izi^i swiy™ cm -3. (12) 

p - 1 8tt m e c 2 7 mi0 \W 



Using our above values for So, B and 7 m> o, we find n e pa 6.1 x 10 4 (r/ro) 2 '° cm 3 . For 
a r pa 0.96 and s — 2, the density profile of the ejecta scales steeply as p oc r~ 27 . 

The associated progenitor mass loss rate is 

M = —n e m p r 2 v w - M yr" 1 (13) 



where we have assumed a nucleon-to-proton ratio of 2 and v w is the velocity of the stellar 
wind. Assuming a typical wind velocity of v w = 1000 km s _1 , we find a progenitor mass loss 
rate of M pa 7.5 x 10~ 6 M yr" 1 for SN2003L. We note that this derived mass loss rate is 
a factor of ~ 30 larger than that of the Type Ic SN 1998bw and SN 2002ap and yet falls at 
the low end of values observed for Galactic Wolf-Rayet (WR) stars (Cappa, Goss & van der 
Hucht 2004), the favored progenitor model for Type Ibc SNe and gamma-ray bursts. 

The total ejecta energy, E, is then given as 



E = _ r 3 -^-° - erg. 14 

r] e ei0 87T V*o/ 

It should be noted that the total energy depends not only on the ratio So, but also on 
assumed value of e ej0 . For e 6i0 = 0.1 (consistent with Li & Chevalier 1999; Berger, Kulkarni 
& Chevalier 2002) we find an ejecta energy of E pa 8.2 x 10 47 (t/t ) ' 89 erg, just a factor of 
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two below that of SN 1998bw and 300 times larger than that of SN2002ap at similar epochs. 
As discussed by Li & Chevalier (1999), a larger # would increase the energy budget further. 

For comparison, we derive the physical parameters implied by Models 2 and 3 here. 
Model 2 shows a slower temporal evolution of radius, r k 5.3 x 10 15 (t/t )°- S5 cm, and 
thus a more rapid decline in velocity, v ps 0.20(t/t )~°' 15 c. The radio bright ejecta scales as 
E 1.0 x 10 48 (t/to) ' 88 erg and the magnetic field evolves with radius as B ps 3.7(r/r )~ ' 99 
G, similar to Model 1. The main difference of Model 2 is that the electron density shows 
a radial dependence of n e ~ 3.7 x 10 4 (r/ro) -1 ' 63 cm~ 3 , significantly shallower than that 
associated with a constant circumstellar wind. The density profile of the ejecta is thus given 
by p oc r~ 10 ' 7 and similar to that found for the Type II SN 1987A (p oc r~ 9 ; Arnett 1988). 
This implies an evolution of the mass loss rate, M ~ 6.8 x 10~ 6 (r/r ) ' 37 M Q yr~ 4 . 

Model 3 shows an even faster deceleration of the ejecta with v ~ 0.20(t/to)~°' 22 c and 
radial evolution of r ps 5.3 x 10 15 (t/to)°' 78 cm. In this case, the energy increases relatively 
slowly as E ps 1.0 x 10 48 (t/t ) ' 34 erg and the magnetic field falls off more steeply as B ps 
3.7(r/r )~ 1M G. The density profile is fixed at s = 2, giving an evolution of n e ps 3.6 x 
10 4 (r/r )™ 2 cm -3 and a constant mass loss rate of M ^ 6.8 x 10~ 6 M Q yr _1 . We note 
that Model 3 predicts a relatively shallow ejecta density profile with p oc r~ 6 - 25 . Figure 6 
compares the physical parameters predicted by Models 1-3. 

6. External Absorption 

In the analysis presented above we did not account for the possible effects of external 
absorption on the observed radio light-curves. For low mass loss rates such as the value we 
derived through SSA modeling of SN2003L (Section 5), it has generally been found that free- 
free absorption (FFA) contributes only minimally to the observed light-curves. Consistent 
with the SSA model, FFA models assume that the intrinsic radio emission spectrum is due 
to non-thermal synchrotron processes. The observable signature of FFA is a steepening of 
the optically thick light-curves at early time (due to contributions from both SSA and FFA 
processes), as was seen in the case of the Type lib SN1993J (Fransson & Bjornsson 1998). 

6.1. Basic FFA Model 

To begin, we fit a basic external absorption model in which the dominant absorption 
process is assumed to be thermal FFA from a uniform ionized circumstellar medium located 
external to the ejecta. Following (Weiler et al. 1986), this FFA dominated model is given 



by 



'"«=Mdk)'(Ty a ' e ^ mJy - (15) 

with 

where a/ and a T describe the temporal evolution of the flux density and optical depth, 
respectively; K\ and K 2 are normalization constants and (5 is the optically thin spectral 
index. 

As shown in Figure 7, the basic FFA model fit (thin solid line) provides a poor match to 
the data, underestimating the slow turn-on of the radio emission with xl ~ 14-0 (two times 
worse than the SSA model fit, § 5). The fitted parameter values are K\ ps 597, (5 ~ —1.0, 
a f « -0.97, K 2 « 1710 and a T w -1.5. 

According to this FFA dominated model, we estimate the predicted mass loss rate of 
the progenitor star using Equation 17 of Weiler, Panagia & Montes (2001). We find a high 
mass loss rate of M m 1.4 x 10~ 2 M yr -1 (for v w = 1000 km s -1 ). This value is three 
orders of magnitude above the typical mass loss rate for Galactic Wolf-Rayet stars, which 
show M k, 10~ 5 M Q yr" 1 with v w 1000 to 2000 km s" 1 (Cappa, Goss & van der Hucht 
2004). 

6.2. Complex FFA Model 

Better light-curve fits are obtained by applying the complex FFA model of Weiler, 
Panagia & Montes (2001) which includes multiple external absorption processes, each char- 
acterized by the distribution of absorbing material (uniform vs clumpy) and the location of 
the material with respect to the radiating electrons (within the CSM vs distant). We note 
that this FFA model also allows for both internal SSA and FFA components to contribute 
to the shape of the observed radio light-curves. 

Applying this model to the SN 2003L radio light-curves, we find that the two dominant 
absorption processes are internal SSA and external FFA due to a clumpy CSM. In the 
notation of Weiler, Panagia & Montes (2001) this corresponds to K 1 m 2.80 x 10 3 , a —1.11, 
(3 ss -1.24, K 3 ^ 7.66 x 10 4 , 6' w -2.15, K 5 ss 1.28 x 10 3 , 5" w -1.60, with all other 
parameters effectively insignificant. As shown in Figure 7, the complex FFA solution (thick 
solid line) provides an reasonable fit to the data with Xr — 4.5 for 96 dof. 
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To estimate the physical parameters of the SN in this model we remove the external 
absorption due to the clumpy CSM component leaving only the intrinsic synchrotron self- 
absorbed spectrum. As shown in Figure 7, at t ~ 2 days the peak flux density is f p ~ 
102.1 mJy at v v tv 22.5 GHz and (3 ~ —1.1. For these values, Equations 2 and 3 give an 
equipartition radius of r cp ps 5.7 x 10 16 cm, an expansion velocity of v ep ~ 14c, and an ejecta 
energy of E cp ps 7.0 x 10 48 erg. With this velocity, the ejecta radius would be r « 2.4 x 10 18 
cm at t ps 65 days, in violation of the VLBA limit of r < 9.1 x 10 17 cm. To match the 
observed VLBA limit, we relax the assumption of equipartition, but this results in a steep 
increase in the total ejecta energy by a factor of ~ 300, giving E ps 3 x 10 51 . Moreover, 
we note that the model is parameterized for the case of non-relativistic ejecta speeds and 
is therefore inconsistent with the implied relativistic velocity. A consistent, non-relativistic 
solution to the Weiler, Panagia & Montes (2001) model could not be found for the case of 
SN2003L. 



7. Modeling the X-ray emission 

X-ray emission in supernovae can be produced by three processes (Fransson, Lundqvist 
& Chevalier 1996): (1) non-thermal synchrotron emission from radiating electrons, (2) ther- 
mal (free-free) bremsstrahlung emission from material in the circumstellar shock and/or the 
ejecta reverse shock, and (3) inverse Compton scattering of photospheric photons by rela- 
tivistic electrons. We examine each of these scenarios in the context of the bright X-ray 
emission of SN2003L. 



7.1. Synchrotron Emission 

Extrapolating the (optically thin) synchrotron emission at t ~ 40 days from the radio 
to the X-ray band as f u oc u 13 with (3 ~ —1.1, we find that the synchrotron emission under- 
predicts the observed X-ray flux by a factor of ~ 50. The discrepancy is significantly larger 
when a synchrotron cooling break, u c , is included, beyond which the spectrum steepens 
by A/3 = —0.5 (see the Appendix). Using Equation A16 together with the magnetic field 
evolution derived in Section 5, the cooling frequency at t ps 40 days is v c ps 1.1 x 10 11 
Hz. Including this synchrotron break and extrapolating from v c to z^x-ray as f u oc z/~ L6 , 
the synchrotron emission in the X-ray grossly under-predicts the observed flux by 5 orders 
of magnitude. We also note that the flat spectrum, (3 ~ —0.5, of the observed X-rays is 
inconsistent with the steep synchrotron spectral index observed in the radio bands, (3 ~ —1.1. 
We therefore conclude that the X-ray emission is not due to synchrotron processes. 
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7.2. Thermal Bremsstrahlung Emission 

In a scenario dominated by thermal bremsstrahlung processes, the X-ray emission is 
produced by the forward shock plowing into circumstellar material and/or the reverse shock 
heating of the ejecta. In this case, the strength of the X-ray emission depends of the density 
of the emitting material. The shocked material then cools by free-free emission processes. 
For an n e oc r~ 2 wind density profile, free-free luminosity from the forward or reverse shock 
is generalized (Chevalier & Fransson 2003; Sutaria et al. 2003) roughly by 



n n34^ ( ^ V / 1000 km/s\ 2 / t \ 1 _ x 

L x ~ 7.5 x 10 34 C L =— - — erg s 1 (17) 

y 10 5 M© yr _1 y \ v w ) \ 40 days J S 1 ; 

where C L is a constant such that C L = C FS = 1 for the forward shock and C L = C RS = 
(n - 3)(n - 4) 2 /(4n - 8) for the reverse shock. For SSA Model 1 (Section 5), we find 
a r ~ 0.96 and thus n ~ 27 and Crs ~ 127. Using the mass loss rate derived from Model 1, 
the combined free- free luminosity from the forward and reverse shocks is therefore Lx,fs+rs ~ 
5.4 x 10 36 erg s _1 , and a factor ~ 10 3 fainter than that observed. On the other hand, if 
we assume the X-ray emission must be dominated by free-free processes, then this implies 
a mass loss rate of 2.7 x 10" 4 M yr -1 (v w = 1000 km s" 1 ), which is inconsistent with 
both the rate predicted by radio SSA models as well as that associated with FFA modeling 
(Section 6). 



7.3. Inverse Compton Cooling 

X-ray emission from supernovae can also be produced by inverse Compton (IC) cooling 
of the relativistic electrons by the optical photons as was claimed to be important for the 
X-ray emission of SN2002ap (Bjornsson & Fransson 2004). In this scenario, the ratio of 
energy densities in magnetic fields, Ub, and photons, C/ p h, is given roughly by 

/radio _ Ub / 1C \ 

(18) 



/x-ray U pil 

where / ra dio = vfv is the optically thin synchrotron flux in the radio band and fx is the 
corresponding X-ray flux. At t 40 days, the observed radio flux was / ra dio ~ 6.7 x 
IQ-ie er g cm -2 s -i anc j X-ray flux was /x-ray ~ 9.2 x 10~ 15 erg cm -2 s _1 , giving a ratio 

Of /radio//x-ray ~ 0.074. 
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Assuming the timescale for Compton cooling is comparable to the epoch Chandra ob- 
servations, tcomp ~ 40 days, and using Equations (14-18) of Bjornsson & Fransson (2004), 
we are able to roughly match the observed X-ray luminosity. Since we see no evidence for IC 
cooling affects within our radio light-curves, we place a lower limit of v\c > 22 GHz. This 
implies constraints on the magnetic field and shock velocity of B > 0.4 G and v > 0.09c, 
consistent with the values predicted from SSA radio analyses (§ 5). The corresponding limit 
on the energy density in magnetic fields is thus Ub > 0.0048 erg cm -3 and is therefore con- 
sistent with the value found through SSA Model 1 (§5). We note that in these estimates, a 
bolometric luminosity of Lboi,42 ~ 15 has been very roughly estimated based on fitting the 
SN2003L optical light-curve (Soderberg et al. 2004). 

While inverse Compton cooling appears to be a feasible process to produce the luminous 
X-ray emission of SN2003L, it predicts a steep spectral index, /3x-ray ~ A-adio ~ — 11> which 
is inconsistent with the observed value of (3 « —0.5 (2 — a). We note, however, that the poor 
photon statistics may be responsible for this discrepancy. 

8. Discussion and Conclusions 

As the first Type Ibc supernova with a radio and X-ray luminosity within an order of 
magnitude of SN1998bw, SN2003L is clearly an unusual event. Here, we summarize our 
findings. Both SSA and FFA models (Sections 5 and 6) show that the energy of the ejecta 
is > 10 48 erg and nearly comparable to that of SN 1998bw at a similar epoch. In Case 1 the 
equipartition and SSA analyses predict a sub-relativistic velocity, v ~ 0.2c, and a relatively 
low mass loss rate, M « 7.5 x 10~ 6 M yr" 1 which are roughly consistent with a scenario 
where inverse Compton scattering produces the observed X-ray emission. In Case 2, the 
complex FFA models reproduce the observed radio light-curves well, yet they imply a highly 
relativistic ejecta speed, v ~ 14c, or a huge energy output of E > 10 51 erg, as well as a 
high mass loss rate of M « 10~ 2 M Q yr _1 which over-predicts the observed X-ray emission 
assuming a thermal bremsstrahlung emission model. 

Our VLBA observation at t ~ 65 days implies an ejecta size which is inconsistent with 
the extreme velocities required by the complex FFA model in Case 2. Only by deviating 
far from equipartition and significantly increasing the total energy of the ejecta can this 
inconsistency be resolved. Moreover, the predicted mass loss rate is at least two orders of 
magnitude above the observed rates from Wolf-Rayet stars. In conclusion, we find Case 2 to 
be infeasible and we therefore adopt Case 1 the most likely scenario. We note that Case 1 
predicts a mass loss rate comparable to the low end of values observed for local Wolf-Rayet 
stars. 
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While similar to SN 1998bw in radio luminosity and ejecta energy, a main distinguishing 
characteristic between SN 2003L and SN 1998bw is their shock velocity: while the radio 
ejecta of SN 1998bw attained mildly-relativistic speeds, that of SN 2003L was clearly sub- 
relativistic as shown in Section 5. This discrepancy is echoed by optical spectroscopy which 
showed photospheric velocities in excess of 3 x 10 4 km s _1 for SN 1998bw (Patat et al. 2001) 
and 10 4 km s" 1 for SN2003L (Valenti et al. 2003; Matheson et al. 2003; Soderberg et al. 
2004). Therefore, we see that for SN 1998bw, mildly-relativistic ejecta was coupled with 
bright radio emission, while for SN2003L the emission is bright despite the modest ejecta 
velocities. The question is thus, why is the radio ejecta of SN 2003L especially luminous and 
energetic? 

One hypothesis considers the distribution of kinetic energy according to ejecta velocity. 
SN 1998bw is a clear case where a significant fraction of the explosion energy has been 
deposited into mildly-relativistic material. For SN2003L, however, the bulk of the energy is 
coupled to sub-relativistic ejecta. Tan, Matzner & McKee (2001) have shown that in some 
cases, high velocity ejecta can be produced from a whip effect racing down the ejecta. The 
observed diversity in the coupling between velocity and energy could therefore be intrinsic to 
the SN Ibc population. A similar situation is observed for gamma-ray bursts which partition 
their energy into ultra-relativistic emission (7-rays) and mildly-relativistic afterglow emission 
(X-rays, optical, radio) by variable fractions. 

Through our ongoing radio surveys of local SNe Ibc and cosmological gamma-ray bursts, 
we continue to map out the energetics of these explosions. Our goal is to understand the 
link between Type Ibc supernovae and GRBs by studying properties of the explosion (ejecta 
velocity and kinetic energy) in addition to those of the progenitor star (CSM density profile 
and evolutionary mass loss rate). Through these further studies, we aim to shed light on the 
illusive bridge between these two classes of cosmic explosions. 

We are indebted to Barry Clark for his generous scheduling of the VLA for this project, 
and to Jim Ulvestad for enabling the VLBA observations. The authors thank Re'em Sari, 
Claes Fransson, Eli Waxman, Dick Sramek and Michael Rupen for useful discussions. Caltech 
SN research is supported by NSF and NASA grants. AMS acknowledges support by the 
National Radio Astronomy Observatory Graduate Summer Student Research Assistantship 
program. Research by EB is supported by NASA through Hubble Fellowship grant HST- 
HF-01171.01 awarded by the Space Telescope Science Institute, which is operated by the 
Association of Universities for Research in Astronomy, Inc., for NASA, under contract NAS 
5-26555. RAC thanks NSF grant AST-0307366. 
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Appendix 



A. The Synchrotron Self- Absorption Model 



Here we describe the SSA model prescription applied in Section 5. The model was 
adapted from the formalism of (Frail, Waxman & Kulkarni 2000), designed to describe 
gamma-ray burst radio emission following the transition to sub-relativistic, adiabatic ex- 
pansion. We replace the assumption of Sedov-Neumann- Taylor (SNT) dynamics (Zeldovich 
& Raizer 2002) with a general parameterization of the shock evolution. This substitution 
enables us to study the early supernova synchrotron emission, while the ejecta is in the 
free-expansion phase. We note that by adopting the scalings given by the SNT dynamics 
into our generalized equations, the formalism of Frail, Waxman & Kulkarni (2000) is fully 
recovered. 



To first order, we assume that the supernova ejecta is undergoing spherical, homologous 
expansion at sub-relativistic velocity, v. At any time, t, the observed synchrotron emission 
originates from a thin shell of radiating electrons with radius, r, and thickness, r/rj (with 
i] ps 10). The electrons are accelerated into a power-law energy distribution, N^) oc 7" p , 
above a minimum Lorentz factor, j m . We adopt the standard assumption that the energy 
density of the ejecta is partitioned between the fraction in relativistic electrons, e e , and 
the fraction in magnetic fields, €3- For simplicity, we denote the ratio of these energies as 
5 = e e /eB and further assume that it evolves as J oc t" 5 . 

The temporal evolution of the radius, velocity, minimum Lorentz factor and magnetic 
field are then parameterized as 



A.l. Physical Assumptions 




(A2) 



(Al) 



(A3) 



(A4) 



where the subscript corresponds to the parameter values at an (arbitrary) reference time, 
t - Here, the indices a r , a 7 , «g and a$ are determined by the hydrodynamic evolution of the 
ejecta, and are constrained according to the assumed evolutionary model (see Section 5.1). 



A. 2. Emission Spectrum 



In characterizing the synchrotron emission spectrum, we use the standard formalism 
of Rybicki & Lightman (1979). Within this framework, the synchrotron power per unit 
frequency emitted by a single electron is given by 



where 7, e and m e are the Lorentz factor, charge and mass of the electron, respectively. The 
critical frequency, z/ crit , is defined as 



(Rybicki & Lightman 1979). Adopting the notation x = (2/3)(z//z/ crit ), the function F{x) 
describes the total synchrotron power spectrum. 



where K 5 / 3 is the modified Bessel function of 2/3 order (Rybicki & Lightman 1979). This 
function peaks at x — 0.29 and decays rapidly for x > 1. 

Applying the temporal scalings for r, (3, B and 7 m (Equations Al - A4), the flux density 
from a uniform shell of radiating electrons is then given by 

fS) = Cf (4 ° r QS)/2 [(1 - e^)]^^ 2 F,{x)F 2 \x) erg/s/Hz/cm 2 (A8) 
where the optical depth, r u (t), is defined 

/ . \ (p-2)a 7 +(3+p/2)a s +a r +a 5 




(A5) 




(A6) 




(A7) 





(A9) 
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and C = [0, 1] parameterizes the sharpness of the spectral break between optically thick and 
thin regimes. We adopt u m = z/ cr i t (7 = 7 m ) as the characteristic synchrotron frequency, 

-J Hz, (A10) 

and for x = (2/3)(v/v m (t)), the functions, F 2 (x) and F 3 (x), are defined 

F 2 (x) = V3 [ X F{y)y^l 2 dy, F 3 (x) = V3 f F{y)y^' 2 dy (All) 
Jo Jo 

where y is a simple integration variable representing the range y = [0,x]. The temporal 

dependencies of F 2 (x) and F 3 (x) can be computed numerically. For i/ C we find the 
following scalings 

F 2 (x) oc t 0A9 ^ +0 - 73 \ F 3 (x) oc t°- 47 ( P -o.i4) (A12) 

and for v 3> neither function evolves. Therefore, using equations A8 through A12, the 
model flux density at any time is strictly determined by three parameters: Cf, C T and z/ m0 : 



47r?7 ' \87r / \m^c 4 7 mi o/ \2irm e c 

^° = 2^< ^ (A15) 

where d is the distance to the source and # is the value of $ at t . Equations A13 - A15 
above show that Cf, C T and v m $ are in turn determined by four physical parameters: B , 
r , 7 m , and # . 



A. 3. Temporal and Spectral Evolution 



We use the above formalism to determine the spectral evolution of the synchrotron 
emission for a given ordering of the synchrotron break frequencies, u a and u m . Here, we 
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define the self- absorption frequency, v a , as the frequency at which the optical depth is unity: 

Tv a = T(U = V a ) = I- 

For completeness, we give the temporal and spectral evolution for higher frequencies near 
the synchrotron cooling frequency, u c , defined as the the frequency above which electrons 
cool efficiently This frequency is given by 

187rm e ce 

Uc = 1^ HZ (A16) 

and is typically located between the radio and optical observing bands during the first few 
years after the supernova explosion. 



A. 3.1. Case 1: v a < v m 

In Case 1, the self-absorption frequency is well below the characteristic synchrotron 
frequency, u a <C v m . In this scenario, the spectrum peaks at u p m u m . It can be shown that 
the temporal scalings associated with the peak are then given by, 



v, 



• p ^Z/ m OCt 2 ^ +QS , f Vp « f Vm OC t 3a r +3a B -a 7 +a ? (Al?) 



Here, u a <C u m so the evolution of u a and f Va depend on the temporal scaling of F 2 (x) 
(Equation A12). We find 



V a OC t( 2 (P- 136 ) Q 7 + (P+264)a s +100a r +100Q3)/(p+164) ^ ^ ^(9a r +8a s -5a 7 +a„ a +3a 5 )/3 (A18) 

where, for simplicity, we use a Va to denote the temporal dependence of v a . For all frequencies, 
the temporal and frequency dependence of f v is then generalized by 



' l/ 2j.(9a r +8a B —5a~ f —5a Va +3a s )/3 



V < V , 



a 



fv OC < 



; .l/3j-(9a r +8Q S -5Q 7 +3a 5 )/3 <-,,<- 
z/ -(p-l)/2 t (6a r +(5+p)a s +2(p-2)a 7 +2a 5 )/2 is < V < V 



l/ -p/2^(6a r .+(8-5p)a s +2(p-2)a 7 +2a 5 -4p+2)/2 y < y 

where the temporal scalings associated with the synchrotron cooling frequency are 
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U c OCt- 3aB ~ 2 , f Vc OC t3a r +(4-p)a B +(p-2)a 7 +a a -p+l > (A20) 



Case 2: z/ m < i/ a 

In Case 2, the characteristic synchrotron frequency is below the self-absorption fre- 
quency, u m <C u a . In this case, the spectral peak occurs at u p z/ a . Since F 2 (x) is constant 
for v u m , the scalings simplify to 



U p « I/ a OC t(2(p-2)a 7 +2(3 +P /2)a s+ 2a r+ 2a 5 )/( P +4) ) ^ _ ^ ^ ^-^+5^ )/2 _ (A21) 

In this scenario, the evolution of v m is given by Equation A10 while F Um scales as 

f Um ott 2ar+2aB+ba ~<. (A22) 

For all frequencies, the temporal and frequency dependence of f u is then generalized by 



U oc { 



u l/Z t (Aa r -a B )/2 



I/ -(p-l)/2 t (4a r -a fl +(4+p)a„ )/2 ( ^ A23 ' ) 



z/ -p/2 t (2a r + (l-3p)a s +(2+p/2)a„ a -2p+l) v <v 

where the temporal scalings associated with the synchrotron cooling frequency are 

U c OC t~ 3aB ~ 2 , f Uc OC t (4^ + (3-3 P )a s +(4+ P )a, a -2 P +2)/2 > 
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Table I. VLA radio flux density measurements of SN2003L 



Date 


At 


-P4.9 GHz ± <y 


-^8.5 GHz ± CT 


-^15.0 GHz ± CT 


-^22.5 GHz ± cr 


Array 


(UT) 


(days) 


(/•«Jy) 


(/•«Jy) 


(/•«Jy) 


(A«Jy) 


Config. 


2003 Jan 26.2 


25.2 


— 


743±39 


— 


— 


DnC 


2003 Jan 28.3 


27.3 


— 


810±63 


— 


3179±85 


DnC 


2003 Jan 29.3 


28.3 


— 


848±65 


— 


3119±86 


DnC 


2003 Jan 30.4 


29.4 


— 


966±60 


— 


3122±75 


DnC 


2003 Jan 31.3 


30.3 


— 


1051±62 


— 


3080±86 


DnC 


2003 Feb 1.2 


31.2 


— 


947±52 


— 


3252±70 


DnC 


2003 Feb 2.3 


32.3 


— 


883±53 


— 


2973±99 


DnC 


2003 Feb 3.3 


33.3 


— 


— 


2287±121 


— 


DnC 


2003 Feb 6.3 


36.3 


— 


1147±47 


— 


— 


DnC 


2003 Feb 8.3 


38.3 


— 


1211±39 


— 


— 


D 


2003 Feb 11.6 


41.6 


— 


1483±76 


2686±156 


2956±96 


D 


2003 Feb 14.5 


44.5 


— 


1448±62 


2753±153 


2684±124 


D 


2003 Feb 16.4 


46.4 


— 


1413±55 


2719±144 


— 


D 


2003 Feb 18.4 


48.4 


— 


1546±45 


— 


— 


D 


2003 Feb 22.3 


52.3 


— 


— 


3014±137 


2714±68 


D 


2003 Feb 23.3 


53.3 


— 


1854±41 


— 


— 


D 


2003 Feb 25.3 


55.3 


— 


1969±53 


2799±149 


2838±99 


D 


2003 Feb 26.3 


56.3 


— 


— 


2865±131 


— 


D 


2003 Mar 2.3 


60.3 


— 


2283±51 


2924±132 


— 


D 


2003 Mar 7.3 


65.3 


848±64 a 






2263±73 


D 


2003 Mar 8.4 


66.4 




2351±47 


2857±123 


2486±66 


D 


2003 Mar 10.3 


68.3 




2081±52 


2966±115 




D 


2003 Mar 14.4 


72.4 




2527±52 


2901±150 


2422±70 


D 


2003 Mar 17.4 


75.4 




2673±57 


3055±173 




D 


2003 Mar 20.3 


78.3 




2422±61 




2130±223 


D 


2003 Mar 23.3 


81.3 




2500±43 


2664±118 




D 


2003 Mar 27.4 


85.4 




2776±51 




1833±100 


D 


2003 Apr 1.2 


90.2 




2551±54 






D 


2003 Apr 2.1 


91.1 








1819±141 


D 
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Table 1 — Continued 



Date 


At 

L — 1 V 


1 4.9 GHz - 1 - u 


L 8.5 GHz - 1 - u 


1 15.0 GHz - 1 - u 


1 22.5 GHz - 1 - u 




CUT) 


v uct .y^v 




(u Iv) 


(u Iv) 


(u Tvl 




2003 Apr 10.4 


99.4 





— 


2581±320 


— 


D 


2003 Apr 30.0 


119.0 





2532±75 





1078±173 


D 


2003 May 3.1 


122.1 





2534±80 





1327±114 


D 


2003 May 16.3 


135.3 





2455±96 


1925±282 


1146±165 


D 


2003 May 28.1 


147.0 


1989±41 


2200±47 


1452±235 


941±161 


D 


2003 Jun 4.0 


154.0 


2410±114 


2151±77 








A 


2003 Jun 17.1 


167.1 


2555±83 


2210±68 


1279±226 





A 


2003 Jul 1.1 


181.1 





2262±63 








A 


2003 Jul 9.0 


189.0 


2485±90 


2075±61 


1135±125 





A 


2003 Jul 24.0 


204.0 


2380±81 


2048±66 


1177±166 





A 


2003 Aug 11.9 


222.9 


— 


— 
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Fig. 1. — Optical and radio emission from SN2003L at t ~ 35 days. A Palomar 60-inch R- 
band image (greyscale) shows the SN located within a spiral arm of host galaxy NGC 3506. 
Twenty-five radio (8.5 GHz) map contours are over-plotted from 0.15 mJy (3-sigma) to 4.0 
mJy in equally spaced logarithmic intervals. The radio emission peaks at the optical position 
of SN 2003L and diffuse emission from the host galaxy extends a size of approximately 45 
by 45 arcsec. 
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Fig. 2. — Radio light-curves of for the Type Ic SN2003L were taken with the VLA at 
frequencies 4.9, 8.5, 15.0, and 22.5 GHz between January 2003 and July 2004. The single 
VLBA observation (5 GHz) is shown as an encircled dot. SSA Model 1 (thick solid line), 
and Models 2 and 3 (thin solid line) as described in §5 are over-plotted. 
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Fig. 3. — Spectral indices for the radio emission of SN 2003L. The temporal evolution of 
4.9/8.5 GHz, 8.5/15.0 GHz, and 15.0/22.5 GHz spectral indices are shown and SSA Model 
1 (§5) is over-plotted. We highlight the VLBA observation with an encircled dot. 
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Fig. 4.— VLBA image of SN2003L at 4.9 GHz at t « 65 days after the explosion. We find a 
flux density for the supernova of -F4.9 qh z = 848 ± 64 /iJy within a 3.158 x 1.285 mas beam. 
Contours define 1 a increments of 67 \x Jy. At a distance of 92 Mpc, this unresolved VLBA 
detection places a direct constraint on the size of the expanding ejecta of r < 9.1 x 10 17 cm 
at t ~ 65 days. 
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Fig. 5. — Chandra X-ray image of SN2003L at t « 40 days. The 2-10 keV luminosity of 
SN2003L is Lx ~ 9.2 x 10 39 erg/s/Hz for a power-law spectral model. SN2003L is therefore 
the most luminous X-ray SN Ibc (on this timescale) with the exception of SN 1998bw. 
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Fig. 6. — Physical parameters for SN2003L based on SSA Model 1 (thick solid line), Model 
2 (thin solid line) and Model 3 (thin dashed line) as described in Section 5. Left column: 
the temporal evolution from t — 10 — 800 days is shown for the ejecta energy, shock radius 
and average velocity (top to bottom). Right column: the radial profile of the magnetic field, 
electron number density and mass loss rate are shown from r = 5 x 10 15 to 2.1 x 10 17 cm. 
The scalings appearing at the top of each plot correspond to Model I. 




Fig. 7. — Free-free absorption fits for SN 2003L based on formalism from Weiler et al. (1986) 
(basic model; thin solid line) and Weiler, Panagia & Montes (2001) (complex model; solid 
line). The basic model provides a poor fit to the data while the complex model provides 
a significantly better fit by including additional absorption processes. The intrinsic (de- 
absorbed) synchrotron spectrum of the complex model (dashed line) implies a relativistic 
velocity which violates the VLBA constraint at t ~ 65 days. 



